The width of the rapidity distribution in heavy ion collisions 
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We have studied the widths of the rapidity distributions of particles produced in nucleus-nucleus 
collisions at various center of mass energies and as a function of centrality at SPS energies. We 
show that the width of the rapidity distribution is sensitive to longitudinal flow, velocity of sound 
in the medium, and rescattering of particles. We explore the possibility of distinguishing the initial 
hard scattering regime from final state effects by studying the variation in the width of the rapidity 
distribution of the particles with centrality for various pr values. 
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The particle production in heavy ion collisions is 
studied by measuring the particle density in rapidity 
(Y) or in pseudorapidity (rj). The evolution of the 
pseudorapidity density at mid rapidity with beam 
energy and centrality has been one of the main in- 
terests of study in heavy ion collisions . Its scaling 
with the number of participating nucleons and/or 
with the number of collisions is believed to provide 
information on the dynamics of the particle produc- 
tion 0,0- The pseudorapidity density at mid rapid- 
ity is also related to the entropy density . However 
similar importance has not been given to the width 
of the rapidity distributions of the particles (cry). 
With the advent of large acceptance detectors such 
as in RHIC experiments |4j and the energy scan be- 
ing high on the agenda of the RHIC program, one 
can study the evolution of the width of the rapidity 
distribution as a function of beam energy and cen- 
trality. The width of the rapidity distribution is be- 
lieved to be sensitive to the following physics effects: 
(a) Final state rescattering ||, hence a pr depen- 
dence study of the width may help in estimating the 
value of pt that separates the initial hard scattering 
regime from the later stage in heavy ion collisions 
which is dominated by rescattering. (b) The width 
of the rapidity distribution contains the information 
of longitudinal flow 6]. (c) For a given freeze-out 
temperature, the width of the rapidity distribution 
in the Landau hydrodynamical model is found to be 
sensitive to the velocity of sound in the medium • 

In this brief report, we first make a compilation 
of the existing data on the widths of the rapidity 
distributions of the particles as a function of cen- 
ter of mass energy (^snn) and centrality. Then we 
show how longitudinal flow, velocity of sound and 
rescattering effect the width of rapidity distribution. 
Qualitatively the variation in width of the rapidity 
distribution with center of mass energy and central- 
ity can be understood on the basis of above men- 
tioned processes. 

In Fig. n we have plotted the width of the rapid- 
ity distribution, cry, for ir~ and K as a function 
of ^snn |a H> 111 • The solid curve corresponds to 
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FIG. 1: Variation of the width of the rapidity distri- 
bution for 7T~, K + , and K~ with center of mass en- 
ergy. The symbols are for nucleus-nucleus collisions and 
the dashed curve for pions in nucleon-nucleon collisions. 
The solid curve corresponds to \/ln(^/sNN /2m p ) from 
Ref [Tof. The dashed curve follows the solid curve up to 
^/sNN = 25 GeV. The errors on the data point are small 
and are within the symbol size. 



the theoretical prediction from Ref. [T(j based on 
the Landau model, developed for studying the ra- 
pidity distribution for pions in pp collisions |Tl| . In 
this model the width of the rapidity distribution is 
given as -\/ln( y/sNN/2m p ) . The dashed curve cor- 
responds to the experimentally determined width of 
the rapidity distribution for pions produced in pp 
collisions. Where not provided in the references di- 
rectly, the widths were obtained by fitting Gaussian 
distributions with center at Y — to the rapidity dis- 
tributions. The following observations can be made 
from the figure: (a) The width of the rapidity dis- 
tribution increases with increase in ^/sjvjv for both 
nucleon-nucleon and nucleus-nucleus collisions, (b) 
The width of the rapidity distribution for pp col- 
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FIG. 2: Variation of the width of the pseudorapidity dis- 
tribution of charged particles as a function of percentage 
of cross section for two center of mass energies y^ijvjv = 
17.3 GeV and 8.76 GeV. 



lisions deviates from the predictions based on the 
Landau model for ^snn > 25 GeV. (c) > &k- 
in nucleus- nucleus collisions. This reflects their dif- 
ferent interaction cross sections with other particles 
in the medium, (d) The nucleus-nucleus data for pi- 
ons shows a similar trend as the curve corresponding 
to ^/ln(y / J/v7v/2TO p ) . ov values are comparable to 
ctk+ an d they are higher than ok- values for SPS 
and RHIC energies. 

In Fig.[21we have plotted the widths of the pseudo- 
rapidity distributions of the charged particles {a c h) 
as a function of % cross section for ^snn = 8.76 
and 17.3 GeV We observe that a c h increases as 
we go from central to peripheral collisions. 

To understand the variation in the width of the ra- 
pidity and the pseudorapidity distribution in Fig. ^ 
and Fig. |21 we need to know the various physics pro- 
cesses that effect the width. Here we will qualita- 
tively discuss how the width of the rapidity distribu- 
tion is effected by (a) longitudinal flow, (b) velocity 
of sound in medium, and (c) initial and final state 
rescattering. 

The rapidity distribution can be used to study 
the longitudinal flow [||. In Fig. [3] we have plot- 
ted the rapidity distribution of 7r _ for y/sJjN = 8.76 
GeV |8j- Our calculations from a static isotropic 
thermal emission model where the rapidity density 
is given as 
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FIG. 3: (a)Rapidity distribution for tt~ at y^ijvjv = 8.76 
GeV ||. Comparison to a thermal model calculation 
and a thermal model with longitudinal flow ({/3l) = 0.6) 
(b) Variation of average longitudinal flow velocity, {pi), 
with center of mass energy. 



is shown by the dashed curve. The temperature T is 
taken as 120 MeV [l3j, m is mass of the pion and A 
is the normalization constant. We observe that the 
thermal model fails to explain the width of the ra- 
pidity distribution. After including the longitudinal 
flow within the ambit of Bjorken hydrodynamics as 
discussed in Ref. [f| , m the above thermal model, 
the rapidity distribution of pions is found to be well 
explained (solid curve). The rapidity distribution is 
now given as 
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and the average longitudinal velocity is defined as 
(Pl) = tanh(?y ma:E /2). = 0.6 is found to explain 
the pion data at ^/snn = 8.76 GeV as shown by the 
solid curve in Fig.|3ta). 



3 



140 
120 


- • 8.76 GeV 7t" 

c s 2 = 1/3 
-c s 2=1/5 

C2=1/6 „ . 




100 






80 






60 






40 






20 






n 


I I I 


\ • 



Rapidity (Y) 



FIG. 4: Rapidity distribution of ir~ at center of mass 
energy 8.76 GeV compared to the rapidity spectra ob- 
tained from the Landau hydrodynamical model for <? 3 = 
0.166, 0.2 and 0.33. 
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FIG. 5: Variation in the width of the rapidity distribu- 
tion as a function of impact parameter for ^/sjvjv = 62.4 
GeV for Au on Au collisions from AMPT and HIJING 
models. Also shown width of rapidity distribution from 
HIJING for particles having p T > 2 GeV/c. 



Wc fitted the rapidity distributions for pions from 
y/SNN = 2 to 200 GeV to results from the ther- 
mal model with longitudinal flow to obtain the 
Variation of (/3l) with ^snn is shown in Fig. E£b). 
We find the average longitudinal velocity for pions 
approach a value of 1 at RHIC from a value of 0.3 
at AGS energies. This is indicative of the higher de- 
gree of nuclear transparency attained at RHIC com- 
pared to SPS or AGS. It is observed that the results 
in Fig. |3Ib) show qualitatively a similar trend with 
\J snn as seen for the ay for pions in Fig. ^ This 
indicates that the collective behaviour and the fi- 
nal state interactions of the produced particles in 
nucleus-nucleus collisions play an important role in 
determining the width of the rapidity distribution. 

The width of the rapidity distribution is sensitive 
to the velocity of sound in the medium formed at 
freeze-out 0. Fig. 0] shows the rapidity distribution 
of pions at ^/snn = 8.76 GeV compared to rapidity 
distribution obtained for various values of velocity 
of sound using Landau hydrodynamics. Within the 
ambit of Landau hydrodynamics one can show, with 
certain assumptions Q, that the rapidity distribu- 
tion has the form 

% ~ Const'-^tMl (3) 
dY 

where a = 2w f /(l - c 2 s ), ui f = ln(T f /T ), T f is the 
freeze-out temperature, To is the initial temperature, 
c s is the velocity of sound in the medium. We ob- 
serve that for a T f = 120 MeV 13] and T Q = 230 
MeV (obtained from the total multiplicity ;3]), c 2 = 



1/5 explains the data very well. A c 2 value of 1/6 
overpredicts the data and a c 2 value of 1/3 (ideal 
gas) underpredicts the data. The x 2 values for the 
distributions with c 2 — 1/3, 1/5 and 1/6 are 51, 1, 
7 respectively. This shows that the width of the 
rapidity distribution of data is sensitive to the pa- 
rameter c 2 representing the velocity of sound in the 
medium in the above model. It may be mentioned 
that the results are sensitive to the choice of initial 
and freeze-out temperatures also jjj. 

Our analysis of the rapidity distributions of pions, 
kaons and protons at AGS and SPS energies all re- 
veal the same value of c 2 . ~ 1/5 which explains the 
data. This may indicate some kind of universality of 
the matter formed at the freeze-out stage. It may be 
mentioned that the value of c 2 ~ 1/5 has been found 
to be a characteristic value of the speed of sound for 
a gas of hadrons |14|. c 2 ~ 1/5 indicates that the ex- 
pansion of the system is slower than that in an ideal 
gas scenario (c 2 = 1/3). Thus the system formed in 
heavy ion collisions gets more time to interact and 
to reach thermal equilibrium. 

Sensitivity of the width of the rapidity distribu- 
tion to rescattering effects is studied here by using 
A Multi Phase Transport (AMPT) model |l5]. The 
AMPT model includes both initial partonic and fi- 
nal hadronic interactions. It uses the input parton 
distribution from the HIJING model 0] . We found 
that a c h obtained from the AMPT model increases 
with ^/sjvat. The values at ^/sjvjv = 17.3, 62.4, 130 
and 200 GeV are 1.57, 2.4, 2.6 and 2.8 respectively 
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for the impact parameter range from to 3 fm in 
Au+Au collisions. In Fig.[S]the widths of the rapid- 
ity distributions for charged particles is plotted as a 
function of impact parameter for ^snn = 62.4 GeV 
for Au on Au collisions using AMPT and HIJING 
models. We observe that a c h increases as we go to- 
wards higher values of impact parameter for both 
AMPT and HIJING models. We also observe such 
a trend at lower and higher center of mass energies. 
This trend is qualitatively similar to that observed 
in the data shown in Fig. [21 The variation in the 
width is smaller for AMPT than for HIJING. This 
indicates that final state rescattering has an effect 
on the width of the rapidity distribution. 

Also shown in Fig. [5] is the variation of the a c h 
from HIJING for the particles with px > 2 GeV/c 
as a function of impact parameter. We observe that 
for particles having p T > 2 GeV/c, <r ch decreases as 
we go higher in collision impact parameter. Qualita- 
tively one can think of the following picture, rescat- 
tering leads to more isotropic momentum distribu- 
tions (for example radial flow in a hydrodynamical 
picture) and hence will lead to narrower rapidity dis- 
tributions. Particles with very high transverse mo- 
mentum which are basically coming from the initial 
state will not exhibit such isotropy in the momen- 
tum distribution. The width of their rapidity distri- 
bution is expected to show a different variation with 
ccntrality. pr > 2 GeV/ c was chosen for this study 
as RHIC results on elliptic flow show that hydro- 
dynamical calculations agree with the data for px 
< 2 GeV/c 0. Thus studying a c h as a function 
of centrality for various pt ranges may indicate the 



possibility of finding a value of transverse momen- 
tum at which the initial hard scattering stage can be 
distinguished from the later final state rescattering. 

In summary, we have found that the width of the 
rapidity distribution of particles increases with in- 
crease in yj snn in heavy ion collisions. ctk+ > g k- 
reflects the different interaction cross sections of op- 
positely charged kaons with other particles in the 
medium. It has also been observed in the data that 
the width of the rapidity distribution increases with 
increase in impact parameter. We have shown that 
the width of the rapidity distribution can be effected 
by longitudinal flow, velocity of sound, and rescat- 
tering. Longitudinal flow is found to increase with 
increase in ^snn- Qualitatively the variation of av- 
erage longitudinal flow velocity with y^s^vTV shows 
a similar trend as ay for pions with ^Jsnn- For a 
given freeze-out and initial temperature the width of 
the rapidity distribution is sensitive to the velocity of 
sound in the medium. A value of c 2 ~ 1/5 is able to 
explain the rapidity distribution of particles at AGS 
and SPS energies. This indicates, the expansion of 
the system is slower than compared to in an ideal 
gas scenario (c 2 = 1/3). Thus the system formed in 
heavy ion collisions gets more time to interact and to 
reach thermal equilibrium. The results from a multi- 
phase transport model qualitatively show a similar 
dependence of the width of the rapidity distribution 
on the impact parameter as does the data. We have 
also studied the possible px dependence of the width 
of the rapidity distribution of charged particles as a 
function of impact parameter to get an idea about 
initial hard scattering and final state effects. 
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